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ABSTRACT 


The photoelectron asymmetry parameter 3 in LS-coupling is 
obtained as an expansion into contributions from alternative angular 
momentum transfers j^. The physical significance of this expansion of 
3 is shown to be that: 1) The electric dipole interaction transfers 

to the atom a characteristic single angular momentum j t = ^ 0 * where 
is the photoelectron's initial orbital momentum, whereas 2) angular 
momentum transfers f indicate the presence of anisotropic (i.e., 
term-dependent) interaction of the outgoing photoelectron with the 
residual ion. For open-shell atoms the photoelectron-ion interaction 
is generally anisotropic; photoelectron phase shifts and electric dipole 
matrix elements depend on both the multi pi et term of the residual ion 
and the total orbital momentum of the ion-photoelectron final -state 
channel. Consequently 3 depends on the term levels of the residual ion 
and contains contributions from all allowed values of j^. These findings 
contradict the independent particle model theory for 3, which ignores 
final-state electron-ion interaction and to which our expressions reduce 
in the limiting cases for which only is allowed, namely 1) 

spherically symmetric atoms [e.g., closed-shell atoms] and 2) open-shell 
atoms for which the electron-ion interaction is isotropic [e.g., very 
light elements]. Numerical calculations of the asymmetry parameters and 
partial cross sections for photoionization of atomic sulfur are presented 
to illustrate the theory and to demonstrate the information on electron- 
ion dynamics that can be obtained from the theoretical and experimental 
study of 6 for open- shell atoms. 
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I. INTRODUCTION 

We obtain in this paper explicit expressions, in LS-coupling, 
for the angular distribution of photoelectrons produced by electric 
dipole ionization of an arbitrary open- or closed-shell atom. Our 
treatment is based on the angular momentum transfer expansion for the 
differential photoionization cross section and is intended to pro- 
vide a theoretical framework that allows angular distribution calculations 
comparable in accuracy to the best calculations of total phqtoionization 
cross sections. The formulas we obtain show explicitly the influence 
of anisotropic electron-ion interactions on the electron angular dis- 
tribution and at the same time explain the success of the Cooper-Zare^ 
independent particle model theory in predicting such distributions for 
closed-shell atoms. * For other than the lightest open-shell atoms, 
however, we expect anisotropic electron-ion interactions to produce 
photoelectron angular distributions that deviate significantly from the 
predictions of the Cooper-Zare theory. 

Our conclusions, described above, are contained implicitly in 
the LS-coupling formulas for the angular distribution asymmetry 
parameter $ obtained by lipsky 8,9 and by Jacobs and Burke , 10 whose 
formulas are in principle equivalent to ours. The advantage of the 
angular momentum transfer expansion employed in this paper, however, 
is that such conclusions follow explicitly from our formulation. 
Hartree-Fock calculations of the angular distribution of electrons 
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photoi oni zed from atomic sulfur, a typical open-shell atom, are 
presented to illustrate our theoretical predictions. 

In Section II we summarize the angular momentum transfer formula- 
tion of the differential photoionization cross section. We also exhibit 
how the angular momentum transfer probes anisotropic electron-ion. 
interactions. The formulas in this section depend on the amplitude for 
.photoionization with a particular value of the angular momentum transfer. 
The form of this amplitude in LS-coupling, a main result of this paper, 
is obtained in Section III. In Section IV we illustrate the theory by 
calculating the photoelectron angular distribution of atomic sulfur. 
Lastly, we discuss our conclusions in Section V. A brief report of 
these results has been published elsewhere.^ 
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II. SUMMARY OF THE ANGULAR MOMENTUM TRANSFER FORMULATION 

The ejection of an electron € from an unpolarized atomic target 
A by electric dipole interaction with an incident photon y may be 
represented schematically as 

A ( J oV + Y(j Y =1 > V‘ 1) ^ A+(J C lr C ) + 0) 

The differential cross section for this process can be separated into 
contributions characterized by alternative values of the angular 
momentum transfer, 

J t = j Y ‘ 1 = J C + S " J 0’ (2) 

provided no measurement is made of either the photoelecron spin or the 

orientation of the residual ion. The vector is the angular momentum 

transferred between the unobserved initial and final angular momenta in 

the reaction, i.e., between the total angular momentum 0 Q of the target 

+ 

A and the combined angular momenta of the residual ion A and the 
photoelectron spin s, which we denote J = J r + s. Allowed values of 

L j U 

are determined by conservation of angular momentum J and parity n 
in Reaction '(1) * 

0 = °0 + j Y = J C + 5 + * ■ (3) 

11 = _7r o = V-n 4 (V 

The general form of the differential cross section for Reaction (1) 
is 12 

. ■&=&[! + BP 2 (COS0)]. . _ (5) 
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Here a is the total cross section, 9 is the angle between the axis of 

linear polarization of the incident light and the direction of the 

outgoing photoelectron, and 3 is the asymmetry parameter. The dynamical 

features of the angular distribution are thus contained in 6» which may 

assume values in the range -1-3-2, corresponding to distributions 
.22 

varying from sin 0 to cos 6. (Though Eq. (5). assumes linearly polarized 

incident light, unpolarized, 13 partially polarized, 14 and elliptically 
1 5 

polarized incident light produce angular distributions that may be 
expressed in terms of 3.) 

The resolution of Eq. (5) into contributions corresponding to 
alternative values of j t requires firstly that one determine the allowed 
values of from Eqs. (2 ) - ( 4) . Secondly, each value of is charac- 
ten zed as being either parity favored or parity unfavored, corres- 
ponding to whether the parity change of the target, tt tt , is equal to 

U 

+(-l) J t or -(-1)^ respectively. The total cross section o and the 
asymmetry parameter 3 may then be expressed in terms of cross sections 
a(j t ) and asymmetry parameters 3{j t ) for a particular value of o t » as 
follows: 3 

e = E o(j t ) (6) 

^t* J cs 

fav unf 

oB = T { l a { J Vfav^ J Vfav “ ? a ^Vunf } ‘ 

J cs J t h 

In Eq. (7) we have summed the favored and unfavored values of 

separately, but as seen from Eq. (11) below, Eq. (7) represents 3 
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as a weighted average of Note that while Eqs. (6) and (7) also 

have sums over J (cf. Ref. (3), p. 1981), we do not indicate the 
dependence of a( j and 3 (j t ) on quantum numbers other than j t until 
Section III of this paper. This dependence is hidden in the scattering 
amplitudes in terms of which a(j t ) and g(j^) are given by: 

°U t >fav " ITT! US + (j t )| 2 + |S_(j t )| 2 ], • (8) 

0 


a ^t^unf " 2J q +1 ! S o^'t>i 


(9) 


e(j t ) 


t'fav 


(j t + 2)|S + (j t )| 2 +(j t -l)|S_(j t )| 2 -3[j t (j t +1)] !5 [S + (j t )S_(jt) + +c.c.] 


(2j t +1)[|S + (j t )| 2 + |S_(j t )|^] 


( 10 ) 


^tW- 1 < 1T > 

In these equations, X is the photon wavelength divided by 2 tt and 
"c.c" denotes "complex conjugate." The parity favored cross sections 
and asymmetry parameters, Eqs. (8) and (10), depend on photoionization 
amplitudes J + (j t ), the "±" sign denoting the value of the photoelectron's 
orbital angular momentum, 1 = j^+1. The parity unfavored partial cross 
sections in Eq. (9) depend on the amplitudes S* (j*.), the "o" denoting 

0 t 

l - j^. That the asymmetry parameter for any parity unfavored value of 
j t is -1 independent of dynamics, as indicated in Eq. (11), is discussed 
in Ref. (2). Finally, the LS-coupling form of S^(j t ) in terms of re- 
duced electric dipole matrix elements is derived in Section III. Before 
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continuing with this formal development, however, we discuss' in the 
rest of this section the physical significance of the angular momentum 
transfer and, in particular, its role as a probe of anisotropic 
electrons-ion interactions. 

The physical significance of angular momentum transfer as a 
direct probe of anisotropic electron-ion interactions is illustrated 
in Fig. 1. In this analysis it is convenient to think of the photo- 
ionization process as having two stages, namely, an initial stage A 
of photoabsorption proper, and a subsequent stage B of escape of the 
photoelectron from the rest of the atom. The angular momentum transfer 
is always equal to the difference between the angular momentum input to 
the atom (namely, the angular momentum j = 1 of the electric dipole 
interaction) and the angular momentum output from the atom (namely, the 
photoelectron's final state orbital momentum l). Thus the angular 

■+ + 4 1 

momentum transfer, j. = j - A, is the net angular momentum transferred 

t Y 

to (or deposited in) the target by the photoionization process. (Note 

that since we consider experiments in which the photoelectron's spin 

s is not measured, s is included as part of the angular momentum of the 

residual target.). The allowed values of j^, however, are different in 

the two stages of the photoionization process. 

In the initial stage. A (illustrated in Fig. la) the photoabsorptioh 

imparts j =1 unit of orbital momentum to the photoelectron, which has 
• • • 

initial orbital momentum (in an independent particle model), yielding 
a final orbital momentum, Jl 1 = i Q + Therefore in stage A the 
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angular momentum transferred to the target is 



where the magnitude jj. has the single value jj. = £ Q . Furthermore, 
owing to parity conservation, £' = £ ±1 , and hence jj. = £ Q is a parity 
favored angular momentum transfer. 

During the subsequent escape of the photoelectron in stage B 
additional angular momentum transfers can arise, within the allowed 
range determined by Eq. (2), from anisotropic interactions of the photo- 
electron with the rest of the target. In this report we consider only 
spin-independent interactions in LS-coupling. Therefore the interaction 
in stage B is that between the orbital motion of the photoelectron and 
the net orbital motion of the electrons of the residual ion core, as 
illustrated in Fig. lb. This interaction produces a dynamical coupling 
of the respective orbital momenta V and L . Owing to the resulting 

V 

angular momentum exchanges t between the photoelectron and the core, 

only the total angular momentum L is conserved. (It is because of this 

dependence on L that we call these interactions anisotropic.) In 

-► ■+ 

particular, the photoelectron orbital momentum can change from £' to £ 
during the departure of the photoelectron from the atom, in which case 
the angular momentum transfer is no longer jj. = -£ q but 

“► _v * "V 

J t H jy * 1 = Jt " k. O 3 ) 

->■ ■ 

as illustrated in Fig. lc. Note that even if the magnitudes of V and 

*+• 

L' remain unchanged, a precession (albeit quantized in units of k) abbut 
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the total orbital momentum L Is sufficient to produce a change in the 
magnitude of j^. 

It is at this point that the connection between the present for- 
mulation and that of the Cooper-Zare independent particle model** emerges 
most clearly. The Cooper-Zare model treats the residual ion core as 
a spectator to the photoionization process. That is, stage B is ignored 
altogether, in which case only the single (parity- favored) angular 
momentum transfer J t = i arises. In addition, the amplitudes 
S ± (j t * £ q ) assume limiting forms (cf. Eq. (35) below)which, when 
substituted in Eq. (10), give the Cooper-Zare formula for the asymmetry 
parameter. These points are developed in detail in the following 
sections. 
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III. PHOTO IONIZATION AMPLITUDES S^(J t ) in LS-COUPLING 


The scattering amplitudes S^(j' t ) may be expressed as a sum of 

reduced electric dipole matrix elements, each one corresponding to a 

■+ 3 

given total angular momentum J: 


W 5 U J c s)J cs t l sU t>l c, o J o J Y“ 1) 

= n(?c) 3 

0 



((O c s)J cs £, J-||P C1] 



(14) 


Here, n(S) = 4wc(nw/(3* 2 ) , 3 = (2J+1 5^, a c denotes the set of quantum 
numbers necessary to uniquely specify the initial state, and the minus 
sign indicates that the final state is normalized according to in- 
coming-wave boundary conditions. Our task in this section is to obtain 
the LS-coupling form of the reduced electric dipole matrix element in 
Eq. (14). 

Before specializing to LS-coupling, however, let us consider the 
problem in general. The form of the reduced dipole matrix element in 
Eq. (14) is inconvenient for numerical calculation for two reasons. 
Firstly, the final state ( ( J c s ) I is defined in terms of the 
dissociation channel quantum numbers appropriate to the electron-ion 
system at infinite separation. In general it is more convenient to 
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calculate the electric dipole matrix element for transition to one of 
the electron-ion eigenchannel states (aJ|, where a denotes the eigen- 
channel coupling scheme. Secondly, it is much 

more convenient to calculate real matrix elements, and for this reason 
a transformation to the standing wave representation is desirable. 

For these reasons, we expand the dipole matrix element in 
Eq. (14) as follows:^ 

^ J c s)J cs £ » J ’ll p[1] H a o J o ) = 

x r((J c s)J cs 5,|a) J exp(i6(a)) (aJ| |P^| |a Q d o ) (15) 

a 

Each term in the summation in Eq. (15) comprises three elements: 

(1) The phase factor i’^exp i (a(J jt) + 6(a)) , which effects the change 

V 

from incoming-wave to standing-wave normalization. Here' a(J Jt) 
is the Coulomb phase, 

o(0 c l) = arg r(A+l-1 //F), (16) 

which depends on the binding energy I(J C ) of the residual ion fine 
structure level J c through the photoelectron kinetic energy e 
measured in Rydbergs: 

e = ^ - I(J C ). (17) 

The phase 6(a) is the photoelectron phase shift with respect to 
Coulomb waves in the eigenchannel a and represents the effect of 
short-range electron-ion interactions. . 
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(2) The real transformation coefficients ( (^ c s )^ cs ^ I a )^ which relate 
the eigenchannel coupling scheme to the dissociation channel 
coupling scheme. 

(3) The real, reduced dipole matrix elements (aJ | |P^ ^ | \^ Q ) • 

Thus far Eq. (15) and all preceding equations are exact for electric 
dipole transitions. Approximations must be made, however, in the 
representation of the eigenchannels (a0| and their phase shifts 6(a) 
as well as in the representation of the initial state |a Q J 0 ). We pro- 
ceed in the rest of this section to derive the LS-coupling form of 
Eq. (15) and then to reduce that further by assuming the use of radial 
one-electron wavefunctions appropriate for given term levels of the 
ion core and of the electron-ion system. 

In LS-coupling a Q and a are given by 
a = L S 

000 (18) 
a = (L C A)L (S c s)S. 

Implicit in the definition of a Q is that we have an atomic configuration 

having a single open shell where is the orbital angular momentum 

and N is the occupation number. Similarly a implicitly indicates the 

configuration of the final state after photoionization, which is of the 

form The transformation coefficient in Eq. (15) may be found 

16 17 

either algebraically or graphically to be: 
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((( L c S c )J cs )l l (L c a)L(S c s)S)J = 

(-l) 2J cs + < L c + * +L > + < S c +s+S > J 0 Cs 

0 OS 

L,. s r -0 
c c c 


S J,c S 

cs 



(19) 


Finally, we must evaluate the LS-coupling form of the reduced dipole 
matrix element: 


,C1] 


M|P LIJ IIVo) - (( L c i) L (S c s)S,J||P [1] ||L o S o ,J o ). (20) 

This evaluation may be carried out graphically, 18 but in what follows 
we shall proceed algebraically. 

The first step in the evaluation of Eq. (20) is to make a 
fractional parentage expansion of the initial state, 

Mo*V = _ 2 _ (21) 

L o S o 

Since the ionization process is spin-independent, the second step is to 
split off the geometrical dependence of the matrix element in Eq. (20) 
on spin and total angular momentum quantum numbers: 1 ^ 

- (a^)L(S c s)SJ||pt 1] |!(L 0 yL 0 (S 0 s)S 0 ,J 0 ) 


= (- 1 ) 


L + s +J 0 +1 J J \ J J o 1 

0 1 L L s n 
0 0 


s( s c ,s 0 ) «(s.s } 


( 22 ) 


«L c l)L||P [1] |'|(C olo )L ) 


X 



15 


to 

The third step is'Veduce the matrix element of the electric dipole 

operator to a one-electron matrix element by factoring out the 

20 

geometrical dependence on core and total orbital momenta: 

>D3n fT 


((L c *)L||P^||(L 0 * 0 )L 0 ) 

' V s {-1 ) L c +a o +L+1 L L ( 

x UI|P C 1 ] IIU- 



5 ( L c’ L o> 


(23) 


In this equation the factor is a weight factor due to the presence 

of N equivalent electrons in the initial state. * The last step is 

to factor the reduced one-electron matrix element into its radial and 
22 

angular parts: 


wip [ 1 ] ii* 0 ) = omic [ 1 ] ik)r^ cScL 


where the angular part is 

U || c ™||0 = (- i) a u 


and the radial part is 


el 


l 1 
0 0 0 


(24) 


(25) 


LSI /» 

l cc = 

cl 


dr 'l'f((L c S c )eS,,LS|r)n^ i {n^ 0 ,L 0 S 0 |r) 


(26) 


Note particularly that this radial part is calculated using radial 
wavef unctions dependent dynamically on the angular momentum and spin 
quantum numbers of the initial state, of the final state, and of the 
residual ion core. 


i 
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Putting Eqs. (21)-(26) together, we find for the reduced dipole 
matrix element in Eq. (20): 

((L c J)L(S c s)S.J||P Cl] ||L 0 S o ,J o ) = 


N 5s U n N ' 1 l. r S r 
0 c c 


’*ol>*o L oV ( -^' C 


L„+S +J +1 AA A A AA 

° ° JJLLA 

0 0 0 


(27) 



Finally, substituting Eqs. (15), (19), (20), and (27) into Eq. (i4) and 
noting that (-1) 0 ^ = (-1)^ we find for the scattering amplitude the 


following result: 

W = nWN^^'^cSc.AolH^o) i^exp.1a(J c A) 



■ L c S c L L C S C L A 2 J ® ^ 2 J +j?+l 

x Z exp ifi J L } ( (-ir J cs +s o + * + L 

L 1 


(28) 



Since we have assumed no dependence of the phase shifts and radial dipole 
matrix elements on total angular momentum J, the sum over J in Eq. (28) 
may be performed analytically using the Biedenharn identity 23 to yield 
the desired expression for the scattering amplitude in LS-coupling: 
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W = n(*)N‘ s U;- 1 L c S c ,l 0 |}i; Lo S 0 )r t exp ia(J^) Q(j t .J c .J cs ) 


x M) < L o +L c +1 > + < S c +s+S o> + ( J cs- J o- j t)3 - p 1 *0 

K ' 0 0 \0 0 0 j 


(29) 


. L c S c L L c S c*" A 2 L ° Lc Jt ' 

* *«T l2 i , L 


L l r * A 

0 C 0 


l 1 L J 


where 

■ k A A A A 

Q( i . _J ) = L S J J 

WVJ t* C* CS ; o 0 c 

Though Eq. (29) gives the LS-coupling form of the scattering 
amplitude we see that there is a geometrical dependence on the quantum 

numbers 0 and J relating to the fine structure levels of the ionic 

cs c [ s L 

c c 

core. (We have neglected any dependence of the phase shifts 6 0 

L S L e 

and radial dipole matrix- elements R £ on the fine structure levels.-) 

We consider this dependence on J cs and J c in turn. 

All of the dependence on J in Eq. (29) is contained in the 
o cs -J 0 cs 

geometrical factor (-1) Q(jt ,J c ,J cs^ 1 w ^ 1ch depends additionally 
on quantum numbers that are either fixed for a given ionization process 
(e.g., s,L-,,S ,L ,S ,J ) or enter incoherently in the differential and 
total cross sections (e.g., in addition to O c$ ). The square of this 
factor, with phase +1, enters into the definition of the cross section 
(Eq . (6) ) and the asymmetry parameter (Eq.(7)), each of which involves 
a summation over O c$ . Accordingly it is convenient to define a new 
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quantity, 

O(j t ,0 c ) 2 = l Q(j t .J c ,J cs ) 2 , (31) 

J cs 

which gives the statistical weight with which ionization probability 
for a given is distributed among the possible fine structure 
levels J c , since 


Z Q (j t >V = 1 


(32) 


Note however that in Eq. (29) there is an additional dynamical depen- 
dence on 0 arising from the Coulomb phase a(J.£). Often, though, 
the fine structure separations of the residual ion are not resolved. 
Then cr(J c &) can be taken as independent of J r and the dependence of 
the cross sections and asymmetry parameter on J can be removed 

V 

altogether by application of Eq. (32). 

Having obtained the form of the scattering amplitude in L$- 
coupling in Eq. (29), it is instructive to return to our discussion 
in the last section concerning the role of as a probe of anisotropic 
interactions, as illustrated in Fig. 1. The allowed values of j t are 
those consistent with the triangular relations, (L L j.) and {£,lj t }, 
implied by the first 6j-symbol in Eq. (29). The coupling of the elec- 
tron to the residual ion (cf. Fig. lb) is reflected in the dependence 

L C S C L L r S r L 

of the phase shifts 6^ and dipole matrix elements R £ ^ on the 
total angular momentum of the electron-ion complex. Only when these 
phases and matrix elements do not depend on L (i.e., when the electron- 



19 


ion interaction is isotropic) is restricted to the single value 
j t = J l Q . For in this case the dynamical weight factors in Eq. (29) may 
be extracted from the summation, since 


' L c S c L L c S c L 1 sotr °P 1c 


ex P R cl 


interaction 


7> ex P 15 ei R eJl’ 


(33) 


and the summation over L may be performed analytically: 

/ > r 


l L' 
L 


L o L c **t 




) 


L o L c V 


= K 2 6 (j W- 


(34) 


The scattering amplitude then depends only on the final orbital angular" 


momentum l of the photoelectron* 


i 1 L 


VW * r exp i( “el +6 £ ll“o ( Q 0 „ '"tt 




(35) 


where the proportionality constant indicates that we have not written 
down all the other factors from Eq. (29) which depend on quantum numbers 
that are fixed for a given photoionization process. These 
other factors do not contribute to the asymmetry parameter -in Eq. (10) 
since they occur in both numerator and denominator and thus cancel out. 
Setting in Eq- (10) and substituting the scattering amplitude 

from Eq. (35) leads to the asymmetry parameter 8(0 of the Cooper- 

4 

Zare independent particle model. 



20 


IV. APPLICATION TO SULFUR PHOTOIONIZATION 

To illustrate the theory developed in the last two sections we 
calculate the angular distribution of photoelectrons ionized from 
atomic sulfur according to the reaction, 

S(3p 4 3 P) + y - S + (3p 3 4 S, 2 D, 2 P) + (36) 

For each of the residual ion terms L S we present in Table I the 

allowed values of photoelectron angular momentum £, angular momentum 

transfer j t> reaction parity (where parity change = +1 is favored and 

parity change = -1 is unfavored), and the allowed values of total 

angular and spin momenta for the electron-ion system. We see that 

the S ion term has only the single angular momentum transfer 

2 2 

J t = = 1 but that the D and P ion terms both have other values 

o 

of j t including parity unfavored values. Notice that for the D and 
2 more 

P ion terms the l = 2 states have two or^l lowed values of LS, implying 
that there will be interference between phase shifts belonging to 
different final state channels. 

For conciseness we shall concentrate in what follows on the 
photoionization reaction leading to the D ion term since it shows 
the strongest anisotropic electron-ion interactions. For this ion 
term the scattering amplitudes S^(j t ) in Eq. (29) for the allowed 
values of £ and listed in Table I are: 
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yi) = C le i(as+is(3D)) R £S ( 3 D) (37a) 

5,(1) ■ c e'^V d ( 3 S) * 1 


ed' 


v ed' 
(37b) 


S d (2) = C |V a d[-I e i5 d (3s) R ed ( 3 S) - le i5d(3p) R £d (3P) + 6 d< 3 °> R^dC 3 D) ] 


(37c) 


S d ( 3 ) = - le ifi d( 3 P) Red (3 P ) °)R £d ( 3 D)] 

(37d) 


In Eq. (37) C denotes those constant factors in Eq. (29) that are common 

- o r 

to all channels, s and d denote Z - 0 and Z = 2, and R .( S)> for example 

L S.L 2 

denotes the radial dipole' matrix element for L C S C - D, L = 0, 

and Z = 2. Note that we have ignored the dependence of the Coulomb 
phase shifts a g and on ion core fine structure levels J c , as 
discussed at the end of the previous section. 

2 

The asymmetry parameter for photoionization to the D ion term is 
given by Eqs. (7)-(ll ) : 

3|S d (l)| 2 - 3-2^[S d (1)S s (l) + +c.c.] - 5|S d (2)| 2 + 2|§ d (3)| 2 


6 = 


(38) 


3|s s o)r + 3|s d (Dr + 5|s d (2)r+ 7|s d o) r 

In Eq. (38) the common factor C in Eq. (37) cancels in numerator and 

denominator. As pointed out in the last section, Eq. (38) reduces to 

4 L c S c L 

the Cooper-Zare result for 3 when the phase shifts 6 0 and radial 

L S L 

dipole matrix elements R became independent of L C S C L. It is of 
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interest to see how this occurs for this particular example. Note 

first that S^J'+O and S d (3)-K) in Eqs. (37c) and (37d) when the phase 

shifts and dipole matrix elements become identical. We also 

see that the squared modulus of each of these scattering amplitudes, 

having f A = 1, is non-zero partly because the resulting factors 

cos(6 d ( J S) - <5 d ( J P)), etc. , in the cross terms are not unity. These 

same factors also arise in cross terms of |S d (l)| and are partly 

responsible for changing the value of this modulus from what its 

(non-zero) value would be in a Cooper-Zare model calculation. For 

these reasons we regard the magnitude of phase shift differences 
L S L L $ L‘ 

to k e an ™di cation 0 f the strength of anisotropic 
electron-ion interactions and hence of the validity of the Cooper- 
Zare model for 3.^ 

The scattering amplitudes in Eqs. (37a) and (37b), having ‘ 

J t * = 1, contribute to the cross section and the asymmetry 

parameter 3 whether the phase shifts and matrix elements are identical 
or not. The scattering amplitudes in Eqs. (37c) and (37d), however, 
having j^. t ,& Q = 1, contribute only when the phases and matrix elements 
are different from one another. An index of the strength of angular 
momentum transfers j t f i Q is thus the fraction [a-a(j t =Jl )]/or where 
a is the photoionization cross section and o(j‘ =& ) is the partial 

V U 

cross section corresponding to j t = i . For photoionization to the 
2 

D ion term of sulfur this ratio is expressed as: 
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[c-o(l )]/c 


5[S d (2)[ 2 + 7|S d (3)[ 2 

3|S s (1)| 2 + 3|S d (l) | 2 + 5|S d (2)| 2 + 7|S d (3) | 2 


(39) 


To evaluate Eqs. (37)-(39) we used continuum Hartree-Fock(HF) 

wave functions obtained by solving the equations given by Dalgarno, 

Henry, and Stewart using methods discussed fully by Kennedy and 

Manson. These continuum wave functions depend on both the ionic 

term level and the total orbital angular momentum. Discrete HF 

single-particle orbitals for the neutral atom and for the ion were 

25 

obtained from the tabulation of Clementi. 

For comparison, we have also carried out a Cooper-Zare type of 

26 

calculation employing Herman-S kill man (HS) wave functions. The 
continuum HS wave functions are calculated in the average sulfur 
potential appropriate to the ground configuration as tabulated by 
Herman and Skillman. These wave functions depend neither on the 
ion core level L S r nor on the total angular momentum L and thus the 
phase shifts and radial dipole matrix elements depend only on and 
Eq. (37) reduces to: 


S s O) = C ie i(a s +6s) R es (40a) 
S d (l) = C |V (a d +5d) R £d (40b) 
S d ( 2 ) = S d (3) = 0 


(40c) 
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Notice that since the HS continuum wave functions do not depend on the 
ionic term level the asymmetry parameters for each ion term, when 
plotted versus photoelectron kinetic energy e, are identical. Discrete 
wave functions for both the ion and the atom were taken to be the 
tabulated HS neutral -atom discrete wave functions.' 

ISA 

In Fig. 2 we have plotted HF phase shifts for the ionic 

term level D as a function of photoelectron kinetic energy e. The 
three allowed values of L are listed in Table I. These phase shifts 
differ by as much as 0.7 radian indicating that anisotropic electron- 
ion interactions are significantly large. 

In Fig. 3 we have plotted the three asymmetry parameters corres- 
ponding to the three alternative ionic term levels resulting from 
photoionization of the sulfur atom. Contrary to the Cooper-Zare model, 
these asymmetry parameters are significantly different from one another 
when plotted as a function of photoelectron kinetic energy. In this 
plot the length formula for electric dipole transitions has been used 
since this is the correct one for HF calculations.^ In Table II, 
however, we list calculated HF asymmetry parameters using both length 
and velocity formulas for the dipole matrix elements in order to show 
that for most energies listed the differences between the asymmetry 
parameters for different ion terms are larger than the length and 
velocity difference for a given ion term. V/e also list for comparison 
the $ parameter calculated using HS wave functions and the Cooper-Zare 
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formula for 0 (i.e., Eq. (40)). The HS wavefunctions are quite 
different from the HF wave functions and thus the HS asymmetry param- 
eter does not seem to be an "average 11 of the HF asymmetry parameters 
at Tow energies. 

In Table III we have plotted HF and HS cross sections for 

photoionization of sulfur. Note that the HS cross section is a total 

cross section and would correspond to the sum of the three HF partial 

cross sections at a given photon energy . However, we have plotted the 

HF partial cross sections as functions of photoelectron kinetic 

energy for comparison with Table II. Comparing Tables II and III, we 

see that the largest differences in the asymmetry parameters occur for 

28 

energies 1.5Ry f e . 2 . 1 Ry . This is just before the Cooper minima 

in the cross sections, which occur in the region 2.1 Ry £ c S 2.8Ry. 

The cross sections in the region 1.5Ry S e < 2.1Ry are of the order 
of 10 cm and thus measurement of 0 for the different thresholds 
should be experimentally possible, if not for sulfur then for some 
other element. Simply put, we wish to emphasize that the differences 
we have found between the asymmetry parameters for the different ionic 
term levels are not dependent on being at a cross section minimum. 
Indeed, as seen in Fig.. 3 and Tables II and III there are measurable 
differences between 0{ P+ S) and 0( P-> P) in the energy range 
O.lRy S e $ Q.8Ry, where the cross sections are of the order of 
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Finally in Fig. 4 we examine the influence of angular momentum 

transfers j t ? on the asymmetry parameter and partial cross section 

2 

for photoionization to the D ionic term level. * The solid line 
represents the asymmetry parameter given by Eq. (38) and plotted also 
in Fig. 3. The dashed line represents 6 calculated according to 
Eq. (38) but setting S^( 2 ) = §^(3) - 0. Note that the result is not 
the Cooper-Zare expression for 6 since we still use Eq. (37b) for 
5^(1) and thus the dependence of phase shifts and radial dipole matrix 
elements on L S L is still important. We see from Fig. 4 that values 

V L 

of j t f = 1 lower fj as much as 0.2 units in the neighborhood of 
e 2 l.ORy. The dot-dashed curve is a plot of the ratio [a-o(l)]/a 
given in Eq. (39). Values of f = 1 contribute as much as 8% 
to the partial cross section in the neighborhood of e = 1 . 5Ry. 
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V. DISCUSSION AND CONCLUSIONS 

We have shown that anisotropic electron-ion interactions in 
atomic sulfur lead to measurable differences between photoelectron 
angular distribution asymmetry parameters corresponding to alternative 
ionic term levels. Similar effects are expected for most open-shell 
atoms. A measure of the strength of anisotropic electron-ion inter- 
actions is the difference between phase shifts for alternative final- 
state channels. In atomic sulfur these phase shift differences are as 
large as 0.7 radian. A separate study of atomic oxygen 7 found phase 
shift differences of. only 0.2 radian and asymmetry parameters that 
were nearly identical for each ionic term level. However, atomic oxygen 
and the other second row elements are regarded as exceptions, since 
they are too light to have strong interactions, and atomic sulfur Is 
regarded as more typical of open-shell atoms in general. Our choice 
of atomic sulfur for study was purely a matter of convenience. We know 
of no experimental data on photoelectron angular distributions for an 
open-shell atom. We emphasize, however, that we expect the magnitude 
of the difference between asymmetry parameters and the magnitude of 
the cross sections to be experimentally measurable for many open- 

shell atoms. 

For the particular case of atomic sulfur we have found that 
angular momentum transfers j t M 0 . «hich do not arise in the Cooper- 
Zare model, 4 contribute only a small but nevertheless significant 
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amount to the asymmetry parameters and to the cross sections. We 
simply do not know whether this will hold true for other open-shell 
atoms. The contributions to the asymmetry parameter from angular 
momentum transfers j t = l are, however, quite different from those in 
the Cooper-Zare model, which has only = & contributions, since 

t U 

the phase shift differences are so large in the different final -state 
channels (cf. Eqs. (37b) and (40b)). 

For closed-shell atoms our formulas reduce rigorously to those 
of the Cooper-Zare model. Unfortunately nearly all experimental 
measurements of photoelectron angular distributions < known to us are 
for closed-shell atoms. Considering the importance of photoelectron 
angular distributions to such diverse areas as radiation dosimetry 
(e.g., 5-ray spectrum) ^ and the physics of the upper atmosphere 
(e.g., conjugate point phenomena) as well as to theoretical physics, 
as emphasized in this paper, we feel that experimental data on photo- 
electron angular distributions for open-shell atoms would be most 
valuable. 

Lastly, we point out that our formulas for photoelectron angular 
distributions have been derived for any electron-ion coupling scheme, 
but worked out in detail only for LS-coupling. In general the electron- 
ion interaction is best described in an intermediate coupling scheme, 
particularly in semi -empirical calculations. Nearly all ab initio 
atomic calculations, however, use the LS-coupling scheme and it is for 
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these calculations that our formulas have been worked out most fully. 
While we- have calculated phase-shifts and dipole matrix elements in 
HF approximation, other more accurate procedures (e.g., many-body 
perturbation theory, random-phase approximation, etc.) may be used to 
compute these quantities for use in our formulas for the asymmetry 
parameter. Similarly, while we have ignored fine structure splittings 
of the ionic core, these may easily be included in angular distribution 
calculations using our formulas as discussed at the end of Section III. 
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Tab 1 e 1 : Allowed values for. the ion core term level ( L S c ) » 

photoelectron orbital angular momentum (£,), angular 

Jt 


momentum 

transfer 

(j t ) 

, reaction parity (^ n ir r (-l 

U v. 

and 

total 

orbi ta 1 

and 

spin angular 

momenta (LS) 

the 

reaction S(3p 4 

3 P)+hv+S + (3p 3 4 

S, 2 D, 2 P) + e 

c S c 

Jl 

h 


Parity 

LS 

4 S 

0 

1 

+ 1 

(favored ) 


4 S 

2 

1 

+1 

(favored) 

3 d 

2 d 

0 

1 

+ 1 

(favored) 

3 d 

2 d 

2 

1 

+ 1 

(favored ) 

3 d, 3 p, 3 s 

2 d 

2 

2 

-1 

(unfavored ) 

3 d, 3 p, 3 s 

2 d 

2 

3 

+1 

(favored ) 

3 d , 3 P, 3 S 

2 P 

0 

1 

+ 1 

(favored ) 

3 P 

2 P 

2 

1 

+ 1 

(favored ) 

3 d, 3 p 

2 P 

2 

2 

-1 

(unfavored ) 

3 d, 3 p 
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4 3 

Table II. HF asymmetry parameters for the reactions S(3p P) + hv 

S + (3p 3 4 S, 2 D, 2 P) + e' as a function of photoelectron kinetic 
energy e using dipole length (velocity) formula and comparison 
with HS asymmetry parameter. 


e(R.y) 

3( 3 P 


!Ll 

8( 3 P 

-y 

2 P) 

8 ( 3 P 

-h 

2 P) 

HS 8 

0.00 

0.144 

( 

0.054) 

0.176 

7 

0.248) 

0.274 

( 

0.360) 

0.35 

0.05 

0.551 

( 

0.495) 

0.584 

( 

0.614) 

0.705 

( 

0.739) 

0.74 

0.10 

0.769 

( 

0.734) 

0.805 

( 

0.814) 

0.942 

( 

0.950) 

0.96 

0.20 

1.052 

( 

1.048) 

1.090 

( 

1.082) 

1.248 

( 

1.229) 

1.25 

0.40 

1.365 

( 

1.407) 

1.402 

( 

1.402) 

1.594 

( 

1.564) 

1.61 

0.60 

1.543 

( 

1.619) 

1.574 

( 

1.604) 

1.779 

( 

1.761) 

1.78 

0.80 

1.659 

( 

1.761) 

1.670 

( 

1.735) 

1.851 

( 

1.856) 

1.82 

1.00 

1.734 

( 

1.849) 

1.693 

( 

1.782) 

1.793 

( 

1.820) 

1.68 

1.25 

1.772 

( 

1.860) 

1.581 

( 

1.640) 

1.472 

( 

1.488) 

— 

1.50 

1.714 

( 

1.665) 

1.234 

( 

1.131) 

0.847 

( 

0.775) 

0.10 

1.80 

1.408 

( 

0.944) 

0.497 

( 

0.185) 

0.044 

(- 

■0.078) 

-0.24 

2.10 

0.724 

( 

0.016) 

-0.127 

(- 

-0.270) 

-0.250 

(- 

-0.268) 

-0.02 

2.30 

0.211 

(- 

-0.210) 

-0.223 

t- 

■0.207) 

-0.173 

(- 

-0.136) 

0.18 

2.60 

-0.108 

(■ 

-0.029) 

-0.029 

( 

0.076) 

0.094 

( 

0.155) 

0.44 

2.80 

-0.013 

( 

0.189) 

0.170 

c 

0.273) 

0.281 

( 

0.339) 

0.60 

3.00 

0.177 

( 

0.393) 

0.360 

( 

0.447) 

‘ 0.450 

( 

0.500) 

0.74 

4.00 

0.933 

( 

1.003) 

0.972 

( 

0.989) 

0.996 

( 

1.009) 

1.10 

8.00 

1<,576 

( 

1.526) 

1.513 

( 

1.523) 

1.522 

( 

1.526) 

1.54 

15.00 

1.697 

( 

1.629) 

1.614 

c 

1.636) 

1.622 

( 

1.636) 

1.62 

30.00 

1.582 

( 

1.561) 

1.534 

t 

1.564) 

1.537 

( 

1.564) 

1.54 
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Table III . HF cross sections for the reactions S(3p H J P) + hv ■* 
S + (3p 3 4 S, 3 D, 3 P) + e“ as a function of photoelec- 
tron kinetic energy e using dipole length (velocity) 

formula and comparison with HS cross section. All 

-18 2 

cross sections are in units of 10“ cm. 


e(Ry ) 

a ( 3 P 

- 4 S) 

at 3 

P 2 d) 

^ < 3 

P + * 2 P) 

0 ( HS ) 

0.00 

. 13.82 

(8.92) 

27.43 

(19.23) 

20.87 

(15.40) 

58.00 

0.05 

14. 61 

(9.11) 

28.00 

(19.18) 

20.78 

(15.04) 

58.78 

0.10 

15.03 

(9.10)' 

27.46 

(18.42) 

19.73 

(14.04) 

57.24 

0.20 

15.04 

(8.65) 

24.21 

(15.65) 

16.03 

(11.09) 

48.33 

0.40 

12.85 

(6.85) 

15.42 

( 9.38) 

8.43 

( 5.62) 

24.60 

0.60 

9.47 

(4.75) 

8.70 

( 5.05) 

4.01 

( 2.62) 

11.84 

0.80 

6. 27 

(2.97) 

4. 67 

( 2.61) 

1.92 

( 1.23) 

4.73 

1.00 

3.84 

(1.73) 

2.47 

( 1.34) 

0.9 5; 

( 0.61) 

2.00 

1.25 

1.96 

(0.82) 

1.14 

( 0.59) 

0.44 

( 0.27) 


1.50 

0.97 

(0. 38) 

0.56 

( 0.29) 

0.23 

( 0.15) . 

0. 58 

1.80 

0.42 

(0.17) 

0.29 

( 0.17) 

0.15 

( 0.10) 

0.52 

2.10 

0.21 

(0.11) 

0.20 

( 0.15) 

0.13 

( 0.10) 

0.57 

2.30 

0.16 

(0.10) 

0.19 

( 0.16) 

0.13 

( 0.10) 

0.61 

2.60 

0.14 

(0.12) 

0.20 

( 0.18) 

0.14 

( 0.12) 

0.68 

2.80 

0.14 

(0.13) 

0.21 

( 0.19) 

0.15 

( 0.12) 

0.71 

3.00 

0.15 

(0.14) 

0.22 

( 0.20) 

0.15 

( 0.13) 

0.7 5 

4. 00 

0.22 

(0.17) 

0.27 

( 0.23) 

0.17 

( 0.14) 

0.80 

8.00 

0.17 

(0.12) 

0.18 

( 0.15) 

0.11 

( 0.09) 

0.42 

15.00 

0.05 

(0. 04) 

0.06 

( 0.05) 

0.04 

t 0.03) 

0.15 

30.00 

0.01 

(0.01 ) 

0.01 

( 0.01) 

0.01 

( 0.01) 

0.03 
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FIGURE CAPTIONS 

Fig. 1. Illustration of the origin of multiple angular momentum transfers 
in atomic photoionization reactions. See text for discussion. 

L c S c L 7 

Fig. 2. Hartree-Fock d-wave phase shifts <5 £d for the D sulfur ion 
term versus photoelectron kinetic energy e for alternative allowed values 
of L. Solid line, L = 0 [i.e., the state 3p 3 ( 2 D)ed 3 S]; dashed line, 

L = 1 ( 3 P) ; dot-dashed line, L = 2( 3 D). 

Fig. 3. Asymmetry parameters $( P+L S ) for the photoionization reactions 
S(3p^ 3 P) ■+ S + (3p 3 L C S C ) + e~ as a function of photoelectron kinetic 
energy. Solid line, 4 S ionic term; dashed line, 2 D; dot-dashed, 2 P. 

Fig. 4. Dependence of the asymmetry parameter &( P+ D) and cross section 
3 2 2 

o( P+ D) for the D ion term of sulfur on angular momentum transfers 
j t t as a function of photoelectron kinetic energy. Left-hand scale 
refers to ( 1 ) the solid line denoting 3 and ( 2 ) the dashed line denoting 
3 ( j t = * 0 = 1)» both for the °P->^D transition. Right-hand scale refers 
to the dot-dashed line which denotes the ratio [a-a( j^=l )]/a for the 
3 P+ 2 D transition. 
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